Re-examination of orbital ordering in LaMnOs: 
Electron-electron and electron-lattice interactions 
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The interactions between Cg orbitals in neighboring sites are investigated in LaMnOa by taking 
into account virtual exchange of electrons and phonons. The spin and orbital ordering temperatures 
and the spin wave dispersion relation are calculated. We find that the orbital ordering is mainly 
caused by the electronic interactions and that the Jahn- Teller coupling is much smaller than that 
reported previously. We propose that the elastic constant shows a characteristic change at the Neel 
temperature by the spin and orbital coupling and the higher-order Jahn- Teller coupling. 
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I. INTRODUCTION 

In some classes of transition metal oxides, degeneracy 
of the d orbitals of a transition metal ion remains and 
electrons have a degree of freedom indicating the opcu- 
pied orbital. This is called orbital degree of freedom.El For 
the colossal magnetoresistanre (CMR) observed in per- 
ovskite manganese oxides,Er □ the orbital degree of free- 
dom is considered to play an important role, because the 
gigantic decrease of resistivity is observed in the vicinity 
of the transition from charge and orbital ordered phase 
to ferromagnetic metallic one. A parent compound of the 
CMR manganites, LaMnOa, shows the orbital ordering 
below 780 K associated with the distortion of a MnOe 
octahedron. It has been experimentally confirmed that 
the type of the orbital ordering is of C-typecl where two 
kinds of orbitals are alternately aligned in the xy plane 
and the planes are stacked along the z axis. In addition to 
the orbital ordering, the so-called ^-type antiferromag- 
netic (AF) ordering appears below 145 K, where spins arg 
aligned parallel (antipararell) in the xy plane {z axis) Bli 
It is well recognized that this anisotropi(j-|nfta|gnetic order- 
ing is stabilized by the orbital ordering.u"ll3 

In 3d transition metal compounds with orbital degen- 
eracy, two kinds of mechanisms have been proposed for 
the orbital ordering. One is caused by the superex- 
change (SE)-type interaction between orbitals in differ- 
ent sites. This interaction originates from the virtual 
exchange of electrons_imder the strong on-site electron- 
electron interactionsJlj'Ej Another mechanism of the or- 
bital ordering is based on the cooperative Jahn- Teller 
(JT) effects where the lattice distortion occurs cooper- 
atively and lifts— tbe. orbital degeneracy in the transi- 
tion metal ions.ll3~t3 The effective interaction between 
orbitals in this mechanism is caused by virtual exchange 
of phonons. However, it is usually difhcult to separate 



contributions of these two mechanisms to the orbital or- 
dering. This is because the two mechanisrps, provide the 
effective orbital interactions cooperatively.EHl l3 

It has been supposed that the strong electron-lattice 
interaction exists and is necessary to explain CMR. The 
orbital ordering in LaMaQspwas also studied based on the 
cooperative JT effects. tSllSila The energy splitting of the 
two eg orbitals due to the lattice distortion termed the 
JT energy (Ejt) was estimated to be of the order of leV 
by anaUraiag the orbital ordering temperature, oatical 
spectracJ^Ea and the energy band calculationBjEd We 
note that in these analyses, the electron correlation effect 
was not taken into account properly. Actually, the on-site 
Coulomb interaction between electrons was estimated ex- 
peririftentally to be about 7eV which is much larger than 
Ejt S3 Therefore, it is necessary to re-examine the or- 
bital ordering in LaMnOa by considering both the coop- 
erative JT effect and the SE interaction under the strong 
electron correlation on an equal footing. 

In this paper, we investigate the interactions be- 
tween Cg orbitals in neighboring sites originating from 
the electron-electron and electron-lattice interactions in 
LaMnOa. Magnitudes of these interactions are deter- 
mined through the calculation of the spin and orbital 
ordering temperatures and the spin stiffness constants. 
It is shown rtfe&ti t^jt is much smaller than that in 
the literaturell3'E2rE2l and the orbital ordering is mainly 
caused by the electronic interactions. We find that the 
elastic constant shows a characteristic change at the 
Neel temperature by which the coupling constant of the 
higher-order JT effect is estimated. 

In Sec. II, the model Hamiltonian which describes the 
orbital interactions caused by exchanges of electrons and 
phonons is derived. In Sec. Ill, we introduce the mean 
field approximation in the formulation of the orbital and 
spin ordering temperatures. In Sec. IV, by comparing 
the theoretical results of the ordering temperatures and 
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the spin stiffness constants with the experimental values, 
the magnitudes of the orbital interactions are determined 
numerically. Temperature dependence of the elastic con- 
stants are studied in Sec. V. The last section is devoted 
to the summary and discussion. 

II. MODEL 

We start with the following Hamiltonian which in- 
cludes spin, orbital and lattice degrees of freedom: 

H = He+ He-latt + 'Hlatt + 'Hstr + 'He-str + 'HhiJT ■ (1) 

Tie describes the electronic interactions and consists of 
three terms as 

Tie = Hj + Hh +'Haf- (2) 



Ti.j represents the SE interaction between nearest neigh- 
boring (NN) Cg electrons derived from|-t|he generalized 
Hubbard model with orbital degeneracylid as 




where Ji = tl/{U' ~ I) and J2 = <g/(C/' + / + 2Jh). 
U, U' and / are the intra-, inter-orbital Coulomb inter- 
actions and the exchange interaction for Cg electrons, re- 
spectively, and a relation U ^ U' + I is assumed. Jh 
is the Hund coupling between eg electron and t2g spin 
Si (5* = 3/2), and to is the transfer intensity between 
NN d3z^^r'2 orbitals along the z axis. Energy splitting 
between two eg orbitals due to the JT effect is neglected 
in the denominators of Ji and J2, because this s pl ii tt iaa 
is much smaller than the Coulomb interactions. tSlE^TEHl 
Si is the spin operator of an eg electron with S — 1/2. 
Tj' is defined as r/ — cos( ^m;)Ti^ — sin{^mi)Tix with 
(mxT'my^mz) = (1,-1,0) where I denotes a direction 

of a bond connecting sites i and j. Ti is the pseu- 
dospin operator for the orbital degree of freedom, and 
{Tiz) = +(— )l/2 corresponds to the state where the 
d3z^-r2 (dx^-y^) orbital is occupied by an electron. The 
second and third terms in Eq. (j^) describe the Hund cou- 
pling between eg and t2g spins and the AF SE interac- 
tion (Jaf) between NN t2g spins, respectively. These are 
given by 

Hh + Ti-AF = — Jh ^ Si ■ 5- + Jaf ^ Sj ■ 5j. (4) 

The second and third terms in Eq. describe the 
electron-lattice interaction and the lattice dynamics, re- 
spectively. Here, we consider the displacement of O ions 




FIG. 1. The modes of the distortion of a MnOe octahedron, 
(a) Q2 and (b) Q3 modes. 




FIG. 2. The modes of the bulk strain, (a) Ux and (b) Mz 
modes. 

along the direction connecting NN Mn ions, since the mo- 
tion of O ions along the other directions does not couple 
(■3^1inearly with the eg orbitals. Thus, He-iatt is given by 

'He-latt = —gjT ^ QilTil , (5) 
i l—x,z 

where gjy is the coupling constant. Qi,j. and Qiz are 
the normal modes of the lattice distortion at site i given 

hy Qix — ^ ^i—xx ^" ^iy ^i—yy) ^^d Qiz — 

■^{2viz - 2vi-z z - Vix + Vi^x X - Viy + Vi^y y)^^ V i^ is the 

displacement of an O ion at -t- and a is the lattice 
constant. These normal modes are schematically shown 
in Fig. |l|. The third term in Eq. (|l|), Hiatt, is given by 

n,au = ^J2 ^{pl,Pj:, + <lk<lj:,), (6) 

k ^—x,y.z 

where or, is the normal coordinate of lattice vibration 

Hi; 

with direction of displacement ^ and momentum k, and 
Pj^^ is the canonical conjugate momentum of g^^. qj:^ and 

v^^ satisfy the relation Vi^ ^ ^Y.i e'''''' with 
N being the total number of Mn sites. The frequency 
of the lattice vibration is independent of k and is given 

by ujj: = with ni being the mass of an O ion, since 
only the spring constant K/2 between NN Mn and O 
ions is taken into account. Tistr and Tle-str in Eq. (|l|) 
describe thcr-d-astic energy and electron-strain coupling, 
respectively,!!!! as 

Hs^r^^iul+ul), (7) 
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and 



Tie 



-250 



'^{uxTix + u^Ti^). (8) 



Here, V is the volume of the system and cq is the elas- 
tic constant. The electron-strain coupling constant go 

is related to gjr as go = i\l Vc^9JT- The bulk dis- 
tortions Ux and Uz are represented by the elastic strain 
^Pp'{PiP' = ^tUi^) as Ux — -7^ (syy — Bxx) and Uz = 
-^{2ezz — eyy — Cxx), respectively .113 Schematic pictures 
of the bulk distortions are presented in Fig. || The last 
term of Eq. (|l|) describes the higher-order JT couphng 
given by 

Ti-hiJT = ^-^pj^ ~ Q'ix)'^iz - '^QizQixTix^, (9) 

with coupling constant B. 

Now, we derive the effective Hamiltonian describing 
the inter-pite orbital interaction through the exchange 
of phononllll from He-iatt and Hiatt- 'He-iatt is rewritten 
by using the Fourier transforms of = e 

andT^, = ^E.e-^^^-~'r.zas 



-ikr 



Here, g^^j^ is defined as 
with 



(10) 




(11) 



(12) 



Then, by using the canonical transformation, the linear 
couplings between Cg electrons and lattice distortion are 
eliminated as 

klV 



+ (13) 



where 



1/2 Cx Cy, ^"y) 



2 I \/^^^^ ^y) ^ 3 (^ 



, (14) 



and Cp = cos kpU. is the new phonon coordinate given 
by 



Ik, 



ki 



ki 



* rri 

9hf=l^-kl ' 



(15) 



and pfe^ is the canonical conjugate momentum for (7^^. 
The first and second terms of the right hand side 
in Eq. ( [l3| ) are denoted by Tio-o and Tiiatt, respec- 
tively. Here, we neglect the noncommutability between 
H (Eq. (|l|)) and qt^. Ho-o includes the self-interaction 
of orbital, which does not contribute to the orbital order- 
disorder transition. Therefore, by subtracting this term, 
we obtain the following form: 



with 



A, 



kW 



'^0-0 ~ y^.'^kw'^-ki'^ki' 

kW 



if Cx Cy, ^^'^^ 

2 V 73*^"^^ ~ ^(-c^ - ^y - ^c^) 



(16) 



(17) 



Then, we obtain the effective Hamiltonian for the spin, 
orbital and lattice degrees of freedom in LaMnOs given 
by 



eff 



He 



"Hlatt + Ho-o + He-str + Hgtr + HhiJT- (18) 



III. MEAN FIELD APPROXIMATION 

In order to calculate the orbital ordering temperature 
from Eq. (18), we introduce the mean field approximation 
at finite temperatures. It is experimentally confirmed 
that the orbital order-disorder transition in LaMnOs is of 
the first order but is close to the second order transition; 
a discontittuity of the orbital order parameter at Too is 
negligible.El Therefore, we expect that the higher-order 
JT coupling, which brings about the first order phase 
transition, is much smaller than the linear JT coupling. 
Thus, we neglect HhijT in Eq. ( p^ and calculate Too- 
We will consider TLm.jt in the calculation of the elastic 
constant presented in Sec. V. 

The ^-type AF spia and C-type orbital orderings are 
observed in LaMn03.Em Two sublattices for the orbital 
(spin) ordering are denoted by A and B {a and 6) and 
the following mean fields are introduced: {Sa,bz), {S^ bz)' 
{Ta,Bx) and {Ta,Bz)- The free energy of the system is 
obtained in the mean field approximation as follows: 



Vco 
2 

N 

+2Jl{{TAzf + {TbzV) 

+ J2{2{Saz){Sbz)~{Sazf 
^{{Taz) + {Tbz))} 



[64y{TAx){TBx) + 2jly{TAz){TBz) 

{Sbz?) 



3 



+ y (Jl - 3J2){{Saz)^ + {Sb^)^ + {Saz){Sb,)) 



NjAF{{sLr + {slMsl) + {siy^ 

^ (In + In z'i, + In + In z^) , 



(19) 



where (3 = 1/T. z^^^^ and z^^^^ in the last four terms rep- 
resent the partition functions of spin and orbital given by 
= Trexp(-/3H^(,)) and z\^^^ = Tr exp(-/3H^(5)), 
respectively. T^a(b) mean field Hamiltonian describ- 

ing the spin state in sublattice a (b) as 



By minimizing the free energy JFq with respect to (Tai) 
for I = X, z and (S'az), the following self-consistent equa- 
tions are obtained: 



(Tai) - Tt[Tai exp(-/3?^^) }/4 



(26) 
(27) 



Equations ( [26| ) and ( p7| ) arc numerically solved under 
the conditions of — 2 go \/ N/Vca{TAz) and Ux = 
which are derived from dT/dux — and dT/dux = 0, 
respectively. 



'^a(b) = '^Jxy{Sa(b) z)Sa(b) z + '^Jz{Sb(a) z)Sa(b)z, (20) 

where JJ^^ and Jf represent the effective exchange inter- 
action between NN spins in the xy plane and along the 
z direction, respectively. These are explicitly given by 

Jxy = - 32('^1 ^ "^"^2) 

+ ^{Jl + J2mTAx){TBx) + {Taz){Tbz)) 

- ^J2{{Taz) + {Tbz)) + ^Jaf, (21) 



and 



+ ^UiTAz) + {Tbz)) + ^Jaf- (22) 
16 16 



In Eq. (pO|), we introduce a spin operator 5,; with 5 = 2, 

rewrite Si and S** as Si = jSi and S^* = jSi, respec- 
tively, and eliminate the largest energy parameter Jjj 
in Eq. (p^). Due to the ^-AF spin structure, the rela- 
tion {Saz) — —{Sbz) is satisfied. Ti-^AiB) mean field 
Hamiltonian describing the orbital state in sublattice A 
(B) as 

'^A{B) — ^Jxy{'^B{A)x)TA{B)x + '^JTa(B)z 

+ {'^Jly{TB{A)z) + 4:JI{Ta{B)z))Ta{B)z 

-2go'^^^{uxTA{B)x + UzTa(b)z), (23) 



where J' / n and J are 

xy(z) 



Jly(z) = 7(3 - J2) + (Jl + J2){S.,z){S,+x(^+z)z) (24) 



and 



J ^2J2{{Saz? -{Saz){Sbz)). 



(25) 



respectively. For the observed C-type orbital ordered 
state, we have the following conditions: {Taz) ~ {Tbz) 
and {Tax) = -{Tbx)- 



IV. TRANSITION TEMPERATURES AND SPIN 
WAVE DISPERSION 

Among several parameters in the Hamiltonian Eqs. (|l|) 
and ( |l8| ) , values of Ji , J2 and gjT are determined by cal- 
culating the spin and orbital ordering temperatures and 
the spin wave dispersion relation. The other parame- 
ters are chosen to be Jaf = 1, a^K = 17x 10'' and 
a^co = 2 X 10'' meV, which are derived from the Neel 
temperature in CaMn03,D the phonon fpquency deter- 
mined by thejnfrared absorption spectraE3 and the elas- 
tic constant. LJ The lattice constant a and the static JT 
distortion Q{= ^ Qf^ + Q.^^l^re chosen to be a = 4 A 
and Q = 0.3 A, respectivelyEZl 

Firstly, we calculate the transition temperatures. The 
orbital ordering temperature in the mean field approxi- 
is given by 



mation T^^J 



9jt \ 
K /' 



(28) 



and the Neel temperature for the A-KP ordering Tjy 
given by the solution of the following equation: 



mf 



rpMF 



-8J. 



xy 



4J|, 



(29) 



where J^y and J 



m) and (E2h are the 



functions of T*^^ 



and T^^^ to the ex- 



defined in Eqs. 
By fitting T*^/ 
perimental transition temperatures Too — 780 K (Ref. |^) 
and Tn — 140 K (Refs. | and 0), respectively, J2 and 
9jtQ are calculated as functions of Ji. In general, the 
mean field approximation tends to overestimate the tran- 
sition temperature. Therefore, we also estimate the pa- 
rameter values of Ji , J2 and gjrQ by considering the cor- 
rection of the mean field transition temperatures. For the 
Neel temperature, we revise T^^^ as bT^^^ with b = 0.63, 
which is the ratio between T/v for the 5* = 1/2 AF Hciscn- 
berg model obtained in the high-temperatipe expansion 
and that in the mean field approximation.Ej As for the 
orbital ordering temperature, we revise as aToQ with 
a = 0.75, which is obtained by pthe calculation of T/y 
for the S" = 1/2 AF Ising model.E3 This is because the 
orbital part of the Hamiltonian (Eq. (^)) has a discontin- 
uous symmetry, although this symmetry is higher than 
that of the Ising model. 
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as 2J|y = —3.34 meV and 



The spin stiffness constant provides another conditions 
for Ji, J2 and gjT- Ahhough the JT distortion does not 
directly couple with the spin degree of freedom, qjt mod- 
ifies the orbital state and affects the SE interaction be- 
tween NN spins. We calculate the spin wave dispersion at 
T — Q and compare it with the experimental one. Here, 
the orbital and lattice degrees of freedom are assumed, ta 
be frozen, since the energy scale of orbital excitational3'E£l 
and optical phonon are much larger than that of spin 
wave. Then, the relevant parts of the Hamiltonian in 
Eq. (|l|) are given by 

'HSW —'He+ TLe-latt ■ (30) 

The static distortion of a MnOg octahedron is written as 
{Q^z,Q^.) = g(cos0/^,sin0/^) whcrc = -ei^ = 
The orbital orderedstate is determined in the 
mean field approximation.^ By applying the Holstein- 
Primakoff transformation to the spin operators, the dis- 
persion relation of the spin wave is calculated. Experi- 
mentally, the spin wave in LaMnOs was measured by the 
neutron scattering experiments in Refs. ^ and The 
authors in these papers analyzed the experimental data 
by using the Heisenberg model with the NN SE inter- 
actions. They obtained magnitude of the interaction in 
the xy plane (JJJ and that in the z axis (Jf), which are 
defined in Eqs. (gl]) and ^ 
2 J| = 2.42 meV. 

By fitting the theoretical results of the orbital order- 
ing temperature, the Neel temperature for the A-AF or- 
dering and the spin stiffness constant to the experimen- 
tal values, we estimate Ji, J2 and gjr- J2 and gjrQ 
are plotted as functions of Ji in Fig. || (a) and (b), re- 
spectively. Solid and broken curves are obtained from 
the transition temperatures. Too and Tjv. The actual 
values of the parameters are expected to exist in the 
shaded regions. Dotted curves are obtained from the 
spin stiffness constants. The analyses for Too and T^r 
(solid and broken curves in Fig. || (a) and (b)) show 
that J2 increases and gjxQ decreases with increasing Ji. 
This is because (3Ji — J2) and g^rp/K contribute coop- 
eratively to Too, as shown in Eq. (^). On the other 
hand, the analyses for the spin stiffness constant (dot- 
ted curves in Fig. || (a) and (b)) show that both J2 and 
9jtQ increase with increasing Ji, i.e., Ji competes with 
both J2 and gjxQ- This is attributed to the facts that 
(1) Ji and J2 are the ferromagnetic and antiferromag- 
netic interactions, respectively, and (2) gjxQ favors the 
(d3a;2_r2/d3j,2_^2)-type orbital ordered state where the 
ferromagnetic interaction in the xy plane caused by Ji 
is weaker than that without gjxQ- The magnitudes of 
Ji, J2 and gjxQ are obtained from the intersection points 
of solid (broken) curve and dotted curve in Fig. || (a) and 
(b). We obtain Ji = 75 - 85 meV, J2 = 25 - 40 meV 
and EjT = gjrQ = 50 ~ 100 meV. It is stressed 
that the value of Ejt is much smaller than that in the 
literature, i.e., Ejt > 1 eVMEfEl In the case that 
a^K = 17 X 10'' meV, the coupling constant of the inter- 
site orbital interaction through the exchange of phonons 
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FIG. 3. (a) J2 and (b) qjtQ as functions of Ji. Solid and 
broken curves are obtained from the analyses of transition 
temperatures Too and Tjv. For broken curves, the correction 
from the mean field approximation is taken into account (see 
the text). The actual values of the parameters are expected to 
exist in the shaded regions. Dotted curves are obtained from 
the analysis of the spin stiffness constants. Other parameter 
values are chosen to be Jaf ~ 1 meV, a?K = 17 x 10'' meV, 
B = meV, Q = 0.3 A and a = 4 A. 



g'jrp/K is about 3 10 meV which is much smaller than 
Ji. We conclude that the small Ejt comes from the 
strong Coulomb interaction and the orbital ordering in 
LaMnOa is dominated by the interaction through the vir- 
tual exchange of electrons. If Ji and J2 are neglected 
to estimate Too, we obtai|n_i?jT to be 400 ^ 700 nieV 
for a^K = 17 X lO-^jaeVEl and 600 - 1100 meV for 
a^K = 40 X 10* meV.y The latter value of Ejt is of the 
same order of magnitude given in Rcf. |lj. 



V. ELASTIC CONSTANTS 

The elastic constants jpr pvide information of the 
higher-order JT coupling ,113113 although the coupling con- 
stant B is supposed to be smaller than that in the lin- 
ear JT coupling. In this section, we examine the elastic 
constants taking into account the electron-electron and 
electron-lattice interactions and the higher-order JT cou- 
pling. 

We start with the model Hamiltonian in Eq. (p^). The 
elastic constants are the coefficients of the 5u^ terms in 



5 



the free energy, where Sui is the deviation of strain ui 
from that in the thermal equihbrium. Here, the devia- 
tions of Tai and Tbi from the thermal equilibrium, 5Tai 
and STbi, are also introduced. STm and STbi are induced 
by an external strain 6uii and the relations between them 
are derived later. Now, the free energy is expanded up 
to the second order of STai, STbi and dui as 



5ul) 



N 



QJ^y5TAx5TBx 



2JlySTAzSTBz 



+2JI{{STaz + STbz 



{STa. 



6TBzf} 

ci 



(31) 



Here, J^q is given by Eq. (Rjl) where z\,g-. is replaced by 



-MB) 



Trexp(-/3H^(3) 



A(B) 

)' ^A(B) being given by 



MB) 



MB) 

Vco 



B 



Ngo 



^A(B)z 



A{B)x)Ta{B)z 



—'^Qa{b)zQa(b)xTa(b)x^- (32) 

The last term of this equation comes from the higher- 
order JT coupling. We assume that the equilibrium val- 
ues of Qai appearing in this term are given by Qaz = auz 
and Qax = "^^^{Tax) by considering the definition of 
the strain Uz and the linear JT coupling in Eq. 



CP's 



are the coefficients of T^i in H}^ (Eq. (23)) where {T^ 



I All 



{Tbi) and ui are replaced by {Tai) + 5Tai, {Tbi) + STbi 
and ui + 6ui , respectively. Their explicit forms are given 
by 



GJ^ySTBx 



/^^Oa 
^9o\l —oux, 



Cf = 2Jly5TBz + 'iJlSTAz - 2go^^6uz 
Cf ^ 2JUTaz + ^JISTbz - 2go\I^Suz 



64ySTAx-2go\I^Su 



Vco 



Vco 



N 



(33) 
(34) 
(35) 
(36) 



Here, the term origiiiating from the higher-order JT cou- 
pling is neglected.E3 Fui in Eq. ( pT| ) represents the self- 
correlation function of the orbital given by 



Fw = {Tai){Tai') - Kw 



(37) 



with 



E 

.mm' 



'x{TAl)mm'iTAv)r 



f3 ^ 



(38) 



where represents the m th eigenvalue of H*^ . To derive 
Eq. (13), the conditions {Tax) = -{Tbx). {Taz) = {Tbz), 
Qaz = Qbz and Qax = ~Qbx are used. By using the 
condition dT/dSTai = 0, the following relations between 
Sui's and STaf's are obtained: 



5Tax + 5Tbx 



4/3ffo Vco 



N 



Fxx{l - {Ul - 2Jly)pFzz] + (4J* - 2J*J/3F, 



6Taz - STbz 
STaz + 5Tbz 

5Tax — 5Tbx 
with 

Dx 



4/3go . IVco 
Dx 



N 



■SUxFxz •, 



(39) 
(40) 



Dz \ N 

x[f,,{1 - 64y(3F,x) + 6J*,/3F4}, (41) 



4/3.go Vco 



4/3ffo Vco 



N 



■SuzF^^, 



(42) 



1 - 6JlypFxx){l - {Ul - 24y)pFzz} 
-&Jly{^Jl-2Jly)(3^Fl, 



(43) 



and 



Dz 



(1 + &JlyPFxx)[l - (4Ji + 2Jly)l3Fzz] 

+ Qjly{Ul + 2jly)P^Fl. 



(44) 



By using Eqs. (|39[)-(|42[), the deviation of the free energy 
from the equilibrium value is given by 



T-To^ ]^{cx{T)5ul + Cz{T)5ul). 



(45) 



Note that the term proportional to SuxSuz is absent be- 
cause of the tetragonal symmetry. Cx{T) and Cz{T) are 
the temperature dependent elastic constants for the Ux 
and Uz modes, respectively. Their explicit forms are 
given by 



Cx{T) 

- ^ 
~ Dx 



and 

Cz{T) 

- ^ 
~ Dz 



[l - {6Jly - Ag'o)f3Fxx}{l - (4 J* - 2j',y)(3Fzz} 

(46) 



i6Jly-Agl){4Jl-24y)p'Fl 



(1 + 64yf3Fxx){l ~ (4 J* -I- 2Jly - 4gl)f3Fzz} 



+6Jly{Ul + 2Jly-Agl)P^Fl 



(47) 
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FIG. 4. The elastic constants as functions of T. Solid 
and broken curves denote and Cz, respectively. Param- 
eter values are chosen to be Ji = 75 meV, J2 = 25 meV, 
Jaf = 1 meV, agjT = 6 x 10^ meV, a'^K = 17 x 10'' meV and 
a'' Co — 2 X 10* meV. The value of B is chosen as B = meV 
in (a) and B = 50 meV in (b). 



For T > Too, c.x{T) and Cz{T) are simplified as 

^(3Ji - J2) 



,(T) = c^(T) = Co 



1 9.JT 

2 K 



9l 



i(3Ji- J2) + i^ 



(48) 



The elastic constants are numerically calculated and 
presented as functions of T in Fig. ^. Parameter val- 
ues are chosen to be Ji — 75 meV, J2 = 25 meV, 
Jaf = 1 meV, agjr = 6x10^ meV, a^K = 17x 10'' meV 
and a?co = 2 x lO"* meV. The higher-order JT cou- 
pling B is chosen to be S = meV in Fig. || (a) 
and S = 50 meV in Fig. ^ (b). When B is compara- 
ble to Ji, the discontinuous change in the elastic con- 
stants is found at Too- We also obtain a change in Cx 
at Tm- This change reflects the change of the orbital 
state. As shown in Fig. ||, when B is much smaller 
than Ji (solid curve), the orbital ordered state of the 
■^{dsz^-r^ + dx2^y2/d3z2-r^ - „j^2 )-type above Tn 
changes to the state below T/v where the component of 
the {dy2_z2 /dz^-x^ytype increases. On the other hand, 
when B is comparable to Ji (broken curve), the compo- 



nent of the -^{dsz^- 



3z2- 



- dx2_y2)-tYpe, 
state increases below T^. These changes of the orbital 
states originate from the spin-orbital coupling in Hj in 




defined by 6*^ = tan- {{Tax) / {Taz))- Solid and broken curves 
are the results for B — meV and i? = 50 meV, respectively. 
Other parameter values are the same as those in Fig. ^. The 
insets show the schematic orbital structures in the xy plane. 



Eq. (^). The changes of the elastic constants are under- 
stood as follows: In the -^{d^z^-r^ + dx2-y2 / d^z2-r2 — 
dx2-y2)-type orbital ordered state, {Tax) is almost sat- 
urated as \{Tax)\ = 1/2 and {Taz) = 0. Thus, the ex- 
ternal strain 5ux does not induce 5Tax and Cx is sat- 
urated. The deviation of the orbital state from the 
-^{d'iz2^r2 + dx2^y2ld'iz2^r2 - dx2_y2)-type one causes 
the decrease of Cx as shown in Fig. ^ This fact reminds 
us of the inverse of the parallel spin susceptibility l/xy 
diverging with decreasing T in an antiferromanget. We 
propose that the characteristic change of the elastic con- 
stant at Tat may be used to estimate the coupling con- 
stant B. 



VI. SUMMARY AND DISCUSSION 

We have examined the orbital ordering in LaMnOs and 
the magnitudes of the interactions of the eg orbitals be- 
tween NN sites caused by the virtual exchange of elec- 
trons and phonons. By calculating the orbital and spin 
ordering temperatures and the spin wave dispersion and 
comparing them with the experimental results, we ob- 
tained Ji = 75 85 meV, J2 = 25 40 meV and 
9.jtQ = EjT = 50 '--^ lODptnei/. Ejt is much smaller 
than that in the literaturala'EjEJ which was estimated by 
neglecting the electron-electron interaction. The present 
results indicate that the orbital ordering in LaMnOa is 
mainly caused by virtual exchange of electrons under the 
strong Coulomb interactions. El We calculate the temper- 
ature dependence of the elastic constants by taking into 
account both the electron-electron and electron-lattice 
interactions. It is predicted that the elastic constants 
show the characteristic change at T/v which depends on 
the magnitude of the higher-order JT coupling; when the 
coupling constant B is comparable to Ji (much smaller 
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than Ji), Cx increases (decreases) with decreasing T. 
Through the detailed comparison between theory and ex- 
periment, the value of B may be determined. 

The present results support the recent report on the 
observation of the collective orbital excitation termed or- 
bital wave in LaMnOs by the Raman scattering exper- 
iments. Saitoh et al. have observed three peak, struc- 
tures around 120-160 meV in the Raman spectra.LJ These 
peaks can be explained by neither the two-phonon ex- 
citations nor the magnetic excitations. The theoreti- 
cal results of the Raman spectra from the orbital wave 
agree with the polarization dependence ofj-the spectra 
and their relative intensities in experiment .CJ Since the 
characteristic energy of the orbital excitation is much 
higher than that of the lattice vibration, we introduce 
the adiabatic approximation for the lattice degree of 
freedom in the calculation of the orbital wave. Then, 
the energy of orbital wave is approximately given by 

uJorb = \j (3Ji + ^Ejt){Ji + J2 + ^Ejt). When we 
adopt the parameter values obtained in the present anal- 
yses Ji — 75 meV, J2 — 25 meV and Ejt — 50 meV, 
we obtain LOorb = 187 meV. These numerical values are 
consistent with the observed energies of the Raman shifts 
from the orbital excitations. 
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